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Abstract 

"a^" 

0^ ■ We consider longwave mode of the interface instability in the system comprising of two 

^jqI immiscible fluid layers. The fluids fill out plane horizontal cavity which is subjected 

■ ■ to horizontal harmonic vibration. The analysis is performed within the framework of 



< 



o 

C/3 



"high frequency of the vibration" approximation and the averaging procedure. The 
nonlinear equation (having the form of Newton's second law) for the amplitude of 



^ , interface deformation is obtained by means of multiple scales method. It is shown that 

in ■ 

Q ' (in addition to previously detected quasistationary periodic solutions) the equation has 

o : 

OO ■ a class of quasistationary solitary solutions. 

o 

o\ 



In experimental works by Bezdeneznykh et al. ||T| and by Wolf |^] for a long horizontal reservoir 
filled with two immiscible viscous fluids, an interesting phenomenon was found at the interface: 
the horizontal vibrations lead to the formation of a steady relief. This formation mechanism has 



a threshold nature; it is noteworthy that such a wavy relief appears on the interface only if the 
densities of the two fluids are close enough, i.e. it does not appear for the liquid/gas interface (free 
surface). The interface is absolutely unstable if the heaviest fluid occupies the upper layer; i.e., the 
horizontal vibration does not prevent the evolution of Rayleigh- Taylor instability, in contrast to the 
vertical one which under certain conditions suppresses its evolution. A theoretical description of this 
phenomenon was provided by Lyubimov & Cherepanov ^ within the framework of a high frequency 
(of the vibration) approximation and an averaging procedure; they found that a horizontal vibration 
leads to a quasistationary state i.e., a state where the mean motion is absent but the interface 
oscillates with a small amplitude (of the order of magnitude of the cavity displacement) with 
respect to the steady relief. They also obtained the general equations and boundary conditions for 



^English translation from " Hydrodynamika" , Perm, 1998, pp. 191 - 196 
^e-mail: khenner@math.tau.ac.il 
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Fig. 1. Problem configuration. 



mean and pulsational parts of the fluid velocities. The theory developed in Q made it possible to 
perform the linear stability analysis for the interface. 

In the approach two parameters were assumed to be asymptotically small simultaneously: 
(i) the dimensionless thickness of the viscous skin-layers 5 = h^^^Ju/uj, u being the kinematic 
viscosity and (ii) the dimensionless amplitude of the vibration e = a/h. In this limiting case, the 
possibility of description of parametric resonant effects is absent and only the basic instability mode 
(Kelvin-Helmholtz, of two counter flows) remains. 

The linear stability analysis for inviscid and viscous fluids in the case of finite e and relatively 
low frequencies of the vibration was carried out analytically and numerically in , Q . The trans- 
formation was found which reduces the linear stability problem under inviscid approximation to 
the Mathieu equation. The parametric resonant regions of instability associated with the inten- 
sification of capillary-gravity waves at the interface and the effects due to viscous damping were 
examined. 

In the present work, following the approach of [^, we make the analytical study for longwave 
interface instability in the high-frequency vibrational field. 

1. Let us consider the system of two immiscible, incompressible liquids filling rectangular cavity 
of length L and height h. In the state of rest the heavy liquid (of density pi) occupies the bottom 
region of height hi , and the light liquid (of density P2) - the upper region of height /12 {h = /ii + /i2) . 
We choose Cartesian coordinate system in such a way that the x, y-axis lie in horizontal plane, the 
z-axis is directed vertically, z = corresponds to the unperturbed interface (Fig. 1). Let the cavity 
perform harmonic oscillation along the x-axis, with the amplitude a and frequency u). 

In the basic state (which is a counter flow), the interface could be considered as plane and 
horizontal. For discussion of this issue as well as the approximation of infinite horizontal layer we 
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refer the readers to 0. 

The quasistationary perturbed state is found from the following problem: 



A^' = 0, A^> = 0, 



(1) 



-Hi: ^ = 0; 



Ho 



$ = 0, 



(2) 



z = ^(x) : 



PH2+H1 ^' 



(3) 



(4) 



B 



pH2 + Hi 



pH2 + Hi 



(1+^2)3/2 



const, 



(5) 



p = Pl/p2, 



B 



(?UJ^ ( Pl- ^''^ 
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where ^' and <I> are the streamfunctions of small 2D normal perturbations, a is the coefficient of 
the interface tension, ^ is interface deformation. Equations (l)-(5) are in dimensionless form; the 
length scale is L = [q;/(/5i — P2)5']^^^ and this is also the scale for ^' and ^. The x, z differentiation 
is denoted by the respective subscripts. In the case of equal heights of the layers {Hi = H2 = H), 
the equations (l)-(5) are reduced to equations (2.6)-(2.7) in Q. 



2. Let us consider the amplitude ^(x) of interface deformation to be small, ^ ^ 1. This allows 
to impose boundary conditions on the unperturbed interface, ^ = 0, being accurate up to linear 
terms. According to the main idea of multiple scales method, we introduce the set of lengths, 



xi = ex, 



X2 = e^x, X3 = e^x, 



and we assume that all variables in (l)-(5) are the functions of these lengths. Then we have the 
following expansions: 



df 



dx dxi dx2 



dx^ 



(6) 
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The parameter B, characterizing the vibration intensity, is represented hke B = + e^r, being 
the threshold instabihty value with respect to longwave perturbations, r is super(under)criticality 
parameter. 

From (l)-(5) we get the following relations in the leading order of the expansion in e: 



_ {pH2 + Hif 

2pip-l){Hi + H2y ^""^ 

The equation (8), in case of equal heights of the layers, is reduced to 

_ Hip + lf 

* M.P - 1) ' 

which is the correct instability threshold value 

In the next order we have (instead of ^2 we write just 

^4 = Ciz^ + C2Z^ + C3Z + C4, <^4 = C5Z^ + Cez'^ + C7Z + C8, (9) 



/ Hf + pHl \ 2p{p - l){Hi + H2)\ 3{p - 1) 2 _ c C- const (10) 
^"^^y- 3{pH2 + H,))+ {PH2 + H1Y pH2 + Hi^ C- const. (10) 

Here Ci — Cg are known functions of p, Hi, H2,£,{x). 

It was shown in Q that in thin {H^ < 3) layers of equal thickness H the most dangerous 
are longwave perturbations (in the sense that they appear at the smallest possible destabilization 
amplitude). The same longwave perturbations in the layers of different heights are the most dan- 
gerous if {Hi + PH2) / {pH2 + Hi) < 3. In the following, we consider this case only. Besides, we are 
interested only in the solutions to (10) which are zero at infinity, i.e. for such solutions C = 0. 

The nonlinear equation (10) could be rewritten in the form of Newton's 2nd law: 

^xx = -dU/d^, U=^qre + ^pf, q{p. Hi, H2) > 0, p{p. Hi, H2) > 0. (11) 
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The function U{(,) (potential energy) is presented in Fig. 2, for different values r < (r = 
-0.4, -0.6, -0.75). 

To zero level of U corresponds the deformation of the interface which has the form of quasi- 
stationary soliton (Fig. 3). It's amplitude, as follows from (11), is given by 

The soliton is stable, since the values r < correspond to undercriticality. This gives hope to 
observe it in the future experiments. 

It is noteworthy that besides the solutions of solitary type to eq. (11), in the space of parameters 
we find also periodic solutions, which correspond to values U < 0. These solutions were examined 
in details in 



References 

[1] N.A. Bezdenezhnykh et al., 1984. Control of the fluid interface stability by vibration, electric 
and magnetic fields. /// All-Union Seminar on hydromechanics and heat/mass transfer in 
microgravity. Chernogolovka, p. 18-20 (in russian). 

[2] M.V. Khenner, D.V. Lyubimov, T.S. Belozerova, B. Roux, Stability of plane-parallel vi- 
brational flow in a two-layer system, accepted to European Journal of Fluid Mechanics. 



The preprint, featuring the main results of this work could be found at http://xxx. lanl. gov/ 



( bhysics/99080071 , 04.08.1999) 



[3] D.V. Lyubimov, A. A. Cherepanov, 1987. Development of a steady relief at the interface of 
fluids in a vibrational field. Fluid Dynamics 22, p. 849-854. 

[4] D.V. Lyubimov, M.V. Khenner, M.M.Shotz, 1998. Stability of a fluid interface under tangential 
vibrations. Fluid Dynamics 33(3), p. 318-323. 

[5] G.H.Wolf, 1969. The dynamic stabilization of the Rayleigh -Taylor instability and the corre- 
sponding dynamic equilibrium. Z. Physik 227(3), p. 291-300. 

[6] A.V. Zamarayev, D.V. Lyubimov, A. A. Cherepanov, 1989. On equilibrium interface shapes in 
vibrational field. UB-AS USSR, Hydrodynamics and heat-mass transfer, p. 23-28 (in russian). 



5 



This figure "figure2.gif" is available in "gif" format from: 



http://arXiv.org/ps/patt-sol/9908005vl 



This figure "figure3.gif" is available in "gif" format from: 



http://arXiv.org/ps/patt-sol/9908005vl 



